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1. Introduction 

For a multiple QSO, the propagation time from the source to the observer 
varies from the image i to the image j, and this difference (Ar^) can be 
measured when the source is variable. In general, assuming a flat universe 
without cosmological constant, the parameter Ar^- x Hq (Hq is the Hubble 
constant) depends on the redshifts of the lens and the source, as well as 
the positions of the individual images and the source (0, , 9j , (3) , and the 
scaled surface potential ^f a at Oi and 6j (see, e.g., Blandford and Kundic 
1996, Williams and Schechter 1997). The observations of multiple images 
of the same source are used to infer (3 and the adjustable parameters a = 
(a±,...,a p ) that appear in the picture of the deflector, i.e., a lens model. 
From the lens model corresponding to the lens picture, Ar^-, 9i, 9j and the 
redshifts, one obtains Hq. 

A golden system (which is suitable for determining Hq) must be a mul- 
tiple QSO verifying some properties. The deflector cannot be dark or very 
faint (non-observable), and in order to model the gravitational potential, 
the lensing structure must be simple (e.g., an isolated galaxy). To infer a 
and (3, a relatively large number of constraints is also desirable. On the 
other hand, for obtaining Arjj, some typical scales of intrinsic variability 
of the source should be less than Ar^, and the absence of strong short 
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timescale microlensing is need. The microlensing events mask the intrinsic 
variability. 

The best studied gravitational lens is the double QSO Q0957+561A,B. 
As a result of recent advances in observations and modelling, a ten percent 
measurement of Hq is attainable. We discuss in detail this system (Sect. 2) 
and review the perspectives in a near future (Sect. 3). 

2. Twin QSO (0957+561^,5) 

For Q0957+561A, B, at optical wavelengths, there are epochs of a rapid 
intrinsic variability (a basic condition for measuring the time delay accu- 
rately) and epochs of calmness. They also appears evidences in favour of 
an important short timescale microlensing (Schild 1996), which masks the 
intrinsic variability. So, due to the double behaviour activity /calmness, the 
microlensing events, the gaps in the monitoring (three months per year), 
etc., it is difficult an accurate time delay determination based on a large 
dataset containing several years of observations. 

The rough estimation Atba ~ T = 420 days (Pelt et al. 1996; see also 
the paper by Pijpers in this volume) can be used for making nice datasets. 
A nice dataset contains an active light curve A during a period and 
the light curve B in the interval [£, + T, tf +T]. The microlensing should be 
weak and, on the other hand, the absence of gaps (each summer) is need. 
Very recently, using image A data in the band r (1994 Dec- 1995 May, 
which were kindly provided us by T. Kundic), and our image B data in the 
band R (1996 Feb.-July), we concluded that Atba = 424±3 days (Oscoz et 
al. 1997). The A component was active (Kundic et al. 1995), and also, we 
have not found evidences of strong microlensing in the r — R comparison 
(Oscoz et al. 1997, Goicoechea et al. 1998). A similar result (417±3 days; 
2cr confidence level) has been derived by Kundic et al. (1997). They use 
light curves in the g and r bands. However, a reanalysis based on discrete 
correlation functions and two nice datasets, A r + B r and A g + B g (we 
basically exclude the photometric monitoring of image B in the interval 
1995 Dec-1996 Jan.), shows that 424 days is the most probable delay. For 
every fixed value 6 (days), we construct the function 

N 

5\6) = (l/N)Y / S l [DCC(r l ) - DAC(r t - 6)}\ (1) 
i=i 

where DCC is the discrete A — B cross-correlation function, DAC is the 
discrete A — A autocorrelation function and Si = 1 only when both DCC 
and DAC are defined at Tj and Ti — 6, respectively, and Si = otherwise. 
Then, we search for a minimum #o> such that 6q = Atba- In the two (g — g 
and r — r) comparisons, it is inferred Atba = 424 days (see Fig. 1). 
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Figure 1. DCC (filled circles) vs. DAC (shifted by 424 days; open circles) in the g — g 
comparison. The agreement is excellent (for details, see main text). 



The main deflector is a cD galaxy at the center of a cluster. A good pic- 
ture of the galaxy could be a King-type surface density profile with velocity 
dispersion a v and core angular radius 6 C , plus a point-mass (M^) at the 
centre of the King profile (Falco et al. 1991). Moreover, the gravitational 
effect of the cluster can be represented by means of a quadratic lens with 
a convergence k and a shear 7 with position angle (j). Grogin and Narayan 
(1996) used a = (<j„, C , M bft , 7, 0), P\ = (fax, Ply), Ps = (Psx,P5y) as free 
parameters (Garrett et al. 1994, have fitted the A and B radio images with 
six Gaussian components, being A\ and B\ the respective core components 
and another jet components), since, in spite of the large number of 

constraints, there is a degenerancy in the convergence due the cluster. The 
whole lens model must include the nine parameters inferred from measure- 
ments of the lensed images as well as an estimation of k derived from a 
direct measurement of the mass either in the cD galaxy or the cluster. A 
measurement of a v (light), the ID velocity dispersion of the luminous stars 
in the galaxy, allows to eliminate the cluster degenerancy. 

Falco et al. (1997) derived a v (light) = 279±13 km/s, and so, with the 
whole lens model and the time delay above mentioned, one obtains Hq = 
66^4 km/s/Mpc (2a). This estimate is 10% accurate at la (Oscoz et al. 
1997). We however remark that a new lens model satisfying the constraints 
from radio mapping as well as new optical constraints deduced from HST 
observations is now in progress (Bernstein et al. 1997). Also, the approxi- 
mation a v = a v (light) must be reconsidered (Mediavilla et al. 1998). 
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3. Perspectives in a near future 

The combined effort by several groups of astronomers will lead to the de- 
tailed analysis of a large number of gravitational mirages. For each system, 
if the deflectors are not dark, we must be able to obtain a promising set 
of constraints and some information on suitable pictures of the global lens 
(primary and secondary deflectors). Through a good picture and the con- 
straints, one can infer a good lens model. Moreover, the system must be 
extensively monitored to get the light curves of the different images and 
deduce at least one time delay. This work could be difficult due to proba- 
ble microlensing events, a weak source variability, an unsuitable sampling 
of the image light curves, etc. The technique for determining a time delay 
plays also a role. Even in the old Twin QSO there is an uncertainty (due to 
the methodology) of about one week, which is irrelevant in order to obtain 
Hq from this system, but it could be dramatic in another multiply imaged 
QSO (e.g., the Triple QSO 1115+080Ai - A 2 ,B,C). 

Nowadays, from two individual systems (Twin QSO and Triple QSO), 
it is inferred a mean value of Hq ~ 60 km/s/Mpc, in good agreement 
with other methods (see the paper by C. Frenk in this volume). However, 
new measurements of Hq from gravitational lenses may surprise us. For 
example, when one measures Hq via gravitational lensing, the influence of 
the large-scale structure is not taken into account. The effect introduced 
in the measurement of Hq caused by large clusters and/or large voids, is a 
very interesting topic, which will be soon studied by our group. 
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